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Abstract 

A transient pulse technique is used to 
obtain harmonic forces from a time-marching 
solution of the complete unsteady transonic 
small perturbation potential equation. The 
unsteady pressures and forces acting on a model 
of the NACA 64AQ10 conventional airfoil and the 
MBB A-3 supercri tical airfoil over a range of 
Mach numbers are examined in detail. Flutter 
calculations at constant angle of attack show a 
similar flutter behavior for both airfoils, 
except for a boundary shift in Mach number 
associated with a correspond! ng Mach number 
shift in the unsteady aerodynamic forces. 
Differences in the static aeroelastic twist 
behavior for the two airfoils are significant. 


Nomenclature 

a pitch axis location, referenced to 

midchord, in semi chords 
b semichord length 

CP center-of-pres$ure location 

Cp pressure coefficient 

c’ chord length 

c^ lift coefficient 

C)t h lift coefficient due to plunge 

c ^u lift coefficient due to pitch 

c, n moment coefficient about c/4 

Cin^ moment coefficient due to plunge 

Cm a moment coefficient due to pitch 

g structural damping coefficient 

h plunge displacement, in semichords 

hi plunge amplitude, in semichords 

k reduced frequency, bco/V 

M free stream Mach number 

ri airfoil mass per unit span 

r a radius of gyration, referenced to 

pitch axis, in semi chords 
s Laplace transform variable 

t time 

V free stream velocity 

x streamwise coordinate 

x s shock location 

x a center-of-gravi ty location, refer- 

enced to pitch axis, in semi chords 
z coordinate normal to stream 

»o mean angle of attack, deg 

amplitude of pitch oscillation, deg 
Aa e static elastic twist angle, deg 

ACp pressure coefficient, normalized by 

modal amplitude 
p mass ratio, m/?rpb 2 

P air density 

t nondimensional time, Vt/b 

9 phase angle, deg 

oscillation frequency 
wh uncoupled plunge mode frequency 

uncoupled pitch mode frequency 


Introduction 

The calculation of aeroelastic response in 
the transonic regime continues to be an active 
area of research because of the significance of 
flutter in this speed range. Finite-difference 
codes now allow the analyst to include the non- 
linear aerodynamic effects associated with shock 
wave formation and motion. These codes have 
been most completely developed for the two- 
dimensional, small -perturbation potential flow 
equation. The LTRAN2 code^ and its extensions 
have been widely used for this purpose. 

Edwards, et al.^ summarized applications of 
the finite-difTerence codes to the flutter pro- 
blem. In the present paper, a version of the 
LTRAN2 code which incorporates higher frequency 
effects' 5 is used. 

The aerodynamic forces generated by a time- 
marching, finite-difference code can be incor- 
porated into a flutter analysis in several ways. 
Classical solutions employ the forces for har- 
monic motion in each vibration mode in an eigen- 
value analysis to determine the flutter point. 
The finite-difference codes may be used to gen- 
erate forces for such an application by forcing 
the airfoil to oscillate in the prescribed mode 
at the desired frequency and marching the solu- 
tion in time until the transients have decayed 
and a periodic solution is achieved. This 
method requires that the code be run several 
times to obtain a single flutter point. A 
second approach is to couple the structural and 
aerodynamic equations and time-march the com- 
plete system from an initial disturbance. The 
stability is assessed from the growth or deeav 
of the resulting transient response. Again, 
several runs of the code are required to deter* 
mine the neutral stability point. A third 
approach is to compute the harmonic forces from 
the Fourier transform of the indicial response. 
That is, a step change in the mode is used as 
the initial condition, and the solution is time- 
marched until the transient has decayed. 

A transient pulse technique, 3 which is a 
variation of the indicial response method, is 
used in this paper to obtain the aerodynamic 
forces. The Fourier transform of the resulting 
force response is divided ‘by the transform of 
the airfoil motion to obtain the harmonic trans- 
fer function. A Pade fit of these k-plane 

forces is made to obtain an s-plane representa- 
tion. A root locus solution is then used to 

determine the stability boundary. This proce- 
dure has the advantage that the expensive aero- 
dynamic code need be run only once for each Mach 
number, mean angle-of-attack, and vibration 
mode. Subsequent flutter calculations for vari- 
ations in structural parameters are relatively 
inexpensive. The applicabilty of the linear 
transform techniques to this nonlinear problem 
described above is verified by comparing the 
results with forced harmonic and with transient 
aeroelastic solutions. 


1 





Edwards, £t have shown a striking 

difference in the flutter behavior of a conven- 
tional and a supercritical airfoil in which the 
latter shows a curl back in the flutter speed 
vs. Mach number when static aeroelatic twist is 
included. In this paper, the flutter behavior 
of a model of the NACA 64A010 conventional air- 
foil and the MBB-A3 supercritical airfoil is 
examined in detail for one mean angle of attack. 
The associated unsteady pressures and forces are 
presented. Si mil iarities and differences bet- 
ween the aerodynamics of the two airfoils are 
interpreted to explain the flutter behavior. 


Aerodynamic Method 

All calculations were made with the 
XTRAN2L J time-marching finite-difference code 
which solves the complete unsteady transonic 
small perturbation potential equation. This 
code is an enhanced version of the LTRAN2-NLFr 
code that includes all of the appropriate time- 
dependent terms in the differential equation and 
boundary conditions. The alternating-di rection 
implicit algorithm of Rizzetta and Chin 5 is 
used in a three-time-level scheme to treat the 
second order time derivative in the differential 
equation. The low frequency far-field radiation 
boundary conditions of Engquist and Majda 6 
have been extended to the full frequency equa- 
tion and incorporated into XTRAN2L . 

The XTRAN2L code uses a default 80x61 com- 
putational x-z grid with 51 points on the air- 
foil. The points are uniformly spaced over the 
airfoil with the exception of an additional 
point near the leading edge. The grid extends 
chord lengths upstream and downstream from 
the airfoil and 25 chord lengths above and below 
the airfoil. Reference 3 points out the neces- 
sity of stretching the grid smoothly away from 
the airfoil and of maintaining adequate grid 
resolution in the far field to avoid erroneous 
internal grid reflection of the outgoing waves. 
Details of the XTRAN2L code and the demonstra- 
tion of its accuracy are given in reference 3. 


Results 


Calculations were made for two airfoils at 
Mach numbers from 0.75 to 0.80 in increments of 
0.01. The airfoils were the model of the NACA 
64A010 * ?sted at NASA Ames Research Center 7 

{called 64AMES herein) and the supercri tical 
MBB-A3 of German design. These airfoils are two 
of the AGARD standard configurations. 8 The 
64 AMES airfoil is about 10.6 percent thick and 
has a small amount of camber. The MBB-A3 air- 
foil is an 8.9 percent thick, aft-loaded super- 
critical airfoil. The ordinates of both air- 
foils are given in reference 8. The airfoil 
slopes required by the aerodynamic code were 
obtained from cubic spline fits to the airfoil 
abscissas and ordinates, using an approximation 
to airfoil arc length as the independent 

variable. All calculations were made at mean 
angles of attack of a 0 = 10 f 0 r the 64AMES 
airfoil and a 0 = -0.5° for the M8B-A3 air- 
foil. 


Steady Aerodynamics 

The steady state pressure distributions for 
the two airfoils are shown in Fig. 1. The 
steady angle of attack for each airfoil was 
chosen such that the shock location and lift 
were about equal at the same Mach number. How- 
ever, the shock strength is mu.ch less on the 
aft-loaded supercritical airfoil. Figure 2 
gives the static aerodynamic parameters as func- 
tions of Mach number. As noted above the lift 
coefficients and shock locations are approxi- 
mately equal, although the lift increases more 
rapidly with Mach number on the 64 AMES airfoil. 
As expected, the pitching moments (about the 
quarter chord) and center-of-pressure locations 
are quite different due to the aft loading on 
the supercritical MBB-A3 airfoil. 


Unsteady Aerodynamics 


Pressure distributions and generalized 
aerodynamic forces were computed for harmonic 
oscillations in two modes, pitch about the quar- 
ter chord and vertical translation (plunge). 
All calculations were made for constant mean 
angle of attack. For harmonic oscillation of 
amplitude ai the pitch motion is given by 

a{x ) = a 0 + ql\ sin kx 

where k = bw/V is the reduced frequency based on 
semi chord and t = Vt/b is the nondi mensional 
time. The plunge displacement is given by 

h(i ) = hi sin kx 


where hi is nondimensional plunge amplitude. 

Flutter calculations generally require the 
determi nation of generalized aerodynamic forces 
for a range of frequencies for each structural 
vibration mode. Harmonic airloads may be cal- 
culated with the XTRAN2L code by specifying the 
mode of motion and frequency of oscillation 
Typically four cycles of oscillation are suffi- 
cient for the unwanted transients to die ou + 
Fourier analysis of the last cycle of oscilla- 
tion gives the Fourier harmonic components of 
the response. Each such calculation may require 
as many as 1440 time steps. Because the flutter 
analyst is interested in the forces for a range 
of frequencies, the use of an indicia! method is 
attractive. Here the Fourier transform is used 
to provide the response at all frequencies of 
interest from a single transient response. A 
variation of this approach is used herein. 


ruise technique . - The indicial response is the 
response to a step function change in angle of 
attack. For the equation solved by the XTRAN2L 
code with the second order time derivative 
included, the resulting response includes non- 
physical transients due to the approximation of 
the infinite initial derivative using the 
finite time steps taken. For this reason a 
smoothly varying exponentially shaped pulse 3 
is used. The input pulse in angle of attack is 
given by 


a( t ) = et 0 + (x\e~( t-17.5At)2/4 
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where At is the nondimensional time step (u^ = 
0.5° and at = 5tt/32 were used). 

An example of the input pulse and resulting 
lift and moment response for the 64AMES airfoil 
at M - 0.8 is shown in Fig. 3. The forces have 
not returned to their initial values for the 128 
time steps shown. It is essential that the cal- 
culation be carried out to sufficient time for 
the transient to decay completely for the 
following reason. The harmonic response is 
obtained from the fast Fourier transforms (FFT ) 
of the input and the response. If the final 
value of the time history has not returned to 
the initial value, the low frequency results 
will be seriously in error. In practice 2048 
steps were used, although 1024 steps were some- 
times sufficient. The FFT technique provides 
forces up to the Nyquist frequency of 
k - if /A t = 6.4. The accuracy of the method for 
frequencies as high as k = 2 is demonstrated in 
reference 3. 


Steady results may be obtained from the 
XTRAN2L code using either the steady, successive 
line over-relaxation solver, or the unsteady 
solver with fixed airfoil geometry and time- 
marching to a steady state. These options lead 
to different steady state results as indicated 
in Fig. 4 for the 64AMES airfoil at M = 0.8. 
Here the steady solver was run to convergence 
with a lift coefficient Cj> = 0.3461 obtained. 
The unsteady solver was then used (with the air- 
foil held fixed) which resulted in the transient 
shown in Fig. 4. As was the case for the pulse 
described in the preceding paragraph, it was 
necessary to calculate several thousand time 
steps for this transient to decay. The tran- 
sient shown was for the first 512 time steps. 
The resulting lift coefficient at convergence 
was Ci = 0.3354. It is, therefore, necessary 
to run the unsteady solver with fixed airfoil 
before applying the pulse technique described 
above in order to insure that the initial and 
final conditions for the pulse are the same. 


Harmonic Forces . - The harmonic forces were 
obtained using the pulse technique by performing 
an FFT of the pulse input (Fig. 3(a)) and of the 
force output (Fig. 3(b)) and dividing the latter 
by the former. For a case with 2048 time steps 
this leads to a frequency resolution of 
Ak = 1/160. A typical result is shown in Fig. 
5. The real and imaginary components of the 
four force coefficients are presented for two- 
cases - the 64 AMES airfoil at M = 0.78 and the 
MBB-A3 airfoil at M = 0.79. The most striking 
result is the agreement in the forces for these 
two different airfoils at the different Mach 
numbers. This agreement in the forces with a 
0.01 shift in M holds true across the Mach num- 
ber range for the cases treated in this paper. 
As may be anticipated, this comparison will be 
reflected in the flutter results to be shown 
later. The jagged nature of the curves at low 
frequency is not a serious problem in flutter 
analysis for which k is usually greater than 
0 . 1 . 


obtained from forced harmonic oscillation of the 
airfoil at k = 0.15. The pressure pulse seen 
near 55% chord is a result of the shock wave 
motion. The greatest nonharmonic response 
occurred at the aft end of the shock pulse and 
was as large as the fundamental component for 
a\ = 2°. Away from the shock, the response 
was nearly all at the fundamental, frequency. 

A phase jump of about 120° occurs at the 
aft end of the shock pulse. Away from the 

shock, the phases and normalized magnitudes are 
essentially equal. The width of the pulse 
increases while the height decreases with 
increasing ai in such a way that the overall 
forces remain nearly equal. That is, the force 
coefficients are independent of the amplitude of 
oscillation. The maximum differences in magni- 
tude and phase for the cases shown are about 
0.5% and 1° for lift and 6.4% and 7° for 

moment. 

Harmonic Airloads . - Calculations for 
forced harmonic motion were made for the plunge 
and pitch modes at each Mach number for both 
airfoils at k = 0.15. The plunge amplitude was 
hi = 0.02, and the pitch amplitude was 
ai - 0.5°. The magnitude and phase of the 

lifting pressure coefficients (normalized by 

modal amplitude) for these cases are shown in 

Fig. 7 for the pitch mode and Fig. 8 for the 
plunge mode. The chief difference between the 
two airfoils is in the strength of the shock 
pulses; the location is about the same. This 

result might be anticipated from the correspond- 
ing steady pressures shown in Fig. 1. At the 

same Mach number, the phase jump at the shock is 

somewhat less for the MBB-A3 airfoil than for 

the 64AMES airfoil for the same reason. 

The unsteady forces that result from these 
pressure distributions at k = 0.15 are shown in 
Figs. 9 and 10. Linear theory results from ker- 
nel function calculations for a flat plate air- 
foil are also shown. The lift coefficients 
(Fig. 9) for both airfoils show only mild 
changes with Mach number and agree quite well. 
The moment coefficients do show significant 
changes with Mach number. The most interest!' ng 
observation from these two figures is the agree- 
ment between the two sets of calculations when 
the results for the M8B-A3 airfoil are shifted 
to a 0.01 lower Mach number. That is, the 

forces for the two airfoils are nearly equal 

across the Mach number range if the square sym- 

bols are shifted 0.01 Mach number to the left. 
This conclusion was anticipated in Fig. 5, which 
shows that this result holds over a wide fre- 

quency range. 

In order to verify the accuracy of the 
pulse technique, the forces shown in Figs. 9 and 
10 were compared with those obtained from pulses 
in h and a. The maximum differences in magni- 
tude occurred for the pitch case and were 
0.048 for cp and 0.018 for c m . The 

maximum phase difference was 30 . 0 


Oscillation Amplitude Effects . - The lift- 
ing pressure distribution on the 64AMES airfoil 
at M = 0.78 oscillating in pitch at four differ- 
ent amplitudes is shown in Fig. 6. The first 
harmonic only is shown. These results were 


Flutter 

Flutter boundaries were calculated for each 
airfoil at each Mach number and one mean angle 
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of attack (a 0 = 10 f 0r the 64 AMES airfoil 
and a 0 - -0.5° for the MBB-A3 airfoil). The 
structural parameters were those of Isogai ’s 
Case A 9 which were used in reference 2. The 
elastic-axis location was at a = -2, the center 
of gravity at x a = 1.8, and the radius of 
gyration squared was r a 2 = 3.48. The mass 
ratio was u = 60 and the modal frequencies were 
<*>h = w a = 100 rad/ sec. 

Flutter Solution Method . - The XTRAN2L code 
may be used to obtain aeroelastic solutions by 
time-marching the coupled aerodynamic and struc- 
tural equations. 2 This method requires sever- 
al runs of the code to determine the neutral 
stability (flutter) point. Alternatively, the 
harmonic forces obtained from the pulse techni- 
que described herein may be used in a conven- 
tional V-g flutter analysis. The harmonic 
forces (k -plane) may also be represented in the 
Laplace variable (s-plane), and root locus 
methods may be used to determine the flutter 
behavior. In this paper a Pade fit 1 ^ was used 
to obtain the s-plane representation. It was 
necessary to use care in applying this method 
because of inaccuracies in the k-plane forces at 
very low frequencies (k <0 . 1 ) . There is also 
some question as to the suitability of the 
particular form of fit for transonic airloads. 
In spite of these reservations, the method was 
used successfully by insuring that the Pade fit 
was accurate near the expected flutter fre- 
quency. 

A sample root locus flutter solution for 
the 64 AMES airfoil is shown in Fig. 11. The 
different curves were obtained by varying the 
speed. The symbols represent the time-marching 
aeroelastic solutions. The tic marks on the 

curves occur at speeds which correspond to those 
for the symbols. Excellent agreement is shown, 
which verifies the chain of linearity assump- 
tions used in obtaining the curves, i.e. pulse 
transient to harmonic forces (via FFT ) to root 
locus (via Pade fit). 

Flutter Results . - The flutter boundaries 
for the two airfoils are shown i Fig. 12 as 
flutter speed index V/bu j a /u and flutter 

reduced frequency as functions of Mach number. 
The Mach number range shown covers the region of 
the transonic dip but does not extend past the 
point of minimum flutter speed. The potential 
flow calculations become unreliable beyond about 
M = 0.8 for the 64 AMES airfoil. Difficulties 
with convergence of the aerodynamic code occur 
which appear to correlate with known problems of 
nonuniqueness with potential flow theory. 11 
An example of this problem was observed in a 
calculation for the 64AMES airfoil at M = 0.84 
and a 0 = 00 for which a decidedly nonsym- 
metric pressure distribution was obtained; the 
upper and lower surface shocks were located at 
x/c = 0.82 and 0.69, respectively. For this 
airfoil at M = 0.81 and a 0 = 1° (the angle 
of attack for the flutter calculations) the 
upper surface shock was located at x/c = 0.82 
after 2048 time steps and the solution was not 
converged. 

The flutter boundaries for the two airfoils 
(Fig. 12) are remarkably similar. If the bound- 
ary for the M8B-A3 airfoil is shifted to the 
left by 0.01 Mach number, the curves are nearly 


coincident. This result is to be expected in 
light of the nature of the aerodynamic forces 
shown earlier (Figs. 5, 9, and 10). The simi- 
larity in the flutter boundaries was surprising 
because of the differences in flutter behavior 
for the two airfoils reported in reference 2. 
Additional calculations made to explain these 
differences are described in the * foil owing para- 
graphs. 

Figure 13 shows the present flutter bound- 
aries as lines with the results from reference 2 
and some additional calculations shown as sym- 
bols. The flagged-symbol points, taken from 
reference 2, were calculated using the LTRAN2- 
NLR code 4 and differ from the present XT RAN 2 L 
code results in: (1) equation level (the second 
time derivative term in the governing equation 
was omitted), (2) fini te-difference grid 
(discussed in reference 3), and (3) airfoil 
slope calculation procedure. In addition, the 
NACA 64A010 was used rather than the 64AMES air- 
foil used herein. In an attempt to sort out 
these different effects, several new calcula- 
tions were made. 

The unflagged circles shown in Fig. 13 were 
computed for the NACA 64A010 airfoil using the 
XTRAN2L code. The difference between these 
points and the solid curve for the 64AMES air- 
foil is a thickness effect and is consistent 
with a transonic similarity shift of about 0.007 
in M corresponding to a ratio in thickness of 
1.06 between the airfoils. The difference bet- 
ween the flagged and unflagged circles (both for 
the 64A010 airfoil) is due to the three code 
differences enumerated above. 

The unflagged square symbols for the MBB-A3 
airfoil, shown in Fig. 13, were obtained from 
the XTRAN2L code with the second time derivative 
term set to zero. The rather large difference 
between these symbols and the dashed curve for 
the MBB-A3 airfoil is due solely to the 

difference in equation level (item (1) above). 
The flagged and unflagged squares were computed 
using the same equation level and differ only in 
items (2) and (3) above. This difference is not 
large, but the results from reference 2 
(flagged) do show a steeper transonic dip in 
flutter speed. In particular, at M = 0.8 the 
flutter speed from reference 2 (flagged square) 
has dropped below the speed computed from the 
same equation (unf lagged square) herein. This 
difference in trend emphasizes the sensitivity 
of transonic calculations to the finite- 
difference grid and airfoil slopes used. 

Static Aeroelastic Twist . - The flutter 
results shown in Figs. 12 and 13 were calculated 
holding the mean angle of attack of the section 
fixed. From a strip theory point of view, this 
would require a different root trim angle for 
each point because of the variations in flutter 
speed and in static pitching moment coefficient 
at each point. The elastic twist is computed 
from 

Aa e = V “r “ ~^T ( ) 2 c m ( V M) 

a 

where a 0 is the section angle of attack and 
a r is the unknown root angle required to 
balance the static aeroelastic moment. The 
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pitching moment coefficient c m (a 0> M) depends 
on a 0 ( a 0 = 10 and -0.5<> for the 64 AMES 
and MBB-A3 airfoils, respecti vely ) and on Mach 
number {as shown in Fig. 2 for c m referenced 
to c/4) and is referenced to the twist axis 
here. 

Reference 2 showed a marked difference in 
behavior for the flutter of the conventional and 
supercri tical airfoils with Mach number when the 
root trim angle was held fixed. (See Figs. 15 
and 16 of reference 2 in which the a 0 corre- 
sponds to a r herein). The determination of 
those results required calculation for several 
values of a 0 and interpolation of the result- 
ing root angles to obtain the desired a r . In 
this paper, all calculations were made for con- 
stant airfoil section angle a 0 . Using the 
root locus procedure, it was easy to vary the 
section elastic axis location and the result of 
this variation on the static elastic twist is 
described below. 

In Fig. 14 the elastic twist at flutter is 
shown for each airfoil for three elastic axis 
locations. The value a = -2 was used in all of 
the results shown earlier and the twist angles 
for these curves correspond to the flutter 
boundaries of Fig. 12. Comparing the two air- 
foils for a = -2, the MBB-A3 airfoil has a lar- 
ger twist and shows a greater Mach number effect 
than does the 64AMES airfoil. The aeroelastic 
windup with M of the MBB-A3 airfoil (2.50) is 
about one degree greater than that of the 64AMES 
airfoil (1.5°). This greater twist is a 
result of the larger pitching moment coefficient 
for the supercri tical airfoil (Fig. 2). If cal- 
culations were made for constant a r one would 
expect the flutter si mi 1 iari ties shown in Fig. 
12 to disappear. The shapes of the flutter 
boundaries with Mach number would be quite dif- 
ferent because the airfoil angle of attack a 0 
would change more for the MBB-A3 than for the 
64AMF.S airfoil. 

The value a = -2 places the pitch mode axis 
one semichord length ahead of the leading edge. 
As this axis is moved closer to the leading edge 
(a = -1.8 and -1.6), Fig. 14 shows the resulting 
decrease in elastic twist and a less severe Mach 
number trend. At a * -1.6, the two airfoils 
show nearly the same Mach number effect. For 
this variation the airfoil mass properties were 
held fixed, that is a + x a = -0.2 in all 
cases. 


Conclusions 

1. A transient pulse technique has been shown 
to be an accurate and efficient tool for deter- 
mining aerodynamic forces for use in transonic 
flutter analysis. The necessity of carefully 
monitoring the convergence of the time-marching 
calculation is emphasized. 

2. An s-plane Pade representation for the har- 
monic forces was used to provide a convenient 
and efficient means for performing parameter 
variations in flutter analysis. The technique 
used was sensitive to small inaccuracies in the 
harmonic forces at low frequencies. 


3. Flutter boundaries were calculated for a 
conventional and a supercritical airfoil. The 
mean angles of attack were chosen such that the 
steady lift and shock locations were the same. 
The flutter boundaries were very similar except 
for a shift of about 0.01 in Mach number. 

4. The unsteady forces for the two airfoils 
were compared at a reduced frequency (k=0.15) in 
the flutter range. The lift coefficients were 
quite similar but the moment coefficients showed 
some differences. Again, a shift of 0.01 in 
Mach number brought the forces into reasonable 
agreement. 

5. The unsteady pressure distributions for the 
two airfoils were also compared at a reduced 
frequency of 0.15. Although the shock strengths 
were different, the locations and widths of the 
shock pulses were very similar. 

6. Unsteady pressure distributions for a range 
of pitch oscillation amplitudes showed that 
although the width of the shock pulse increased 
with increased amplitude, the unsteady coeffi- 
cients were nearly constant. 

7. Calculations of the static aeroelastic twist 

at flutter were made for several pitch axis lo- 
cations. The supercritical airfoil showed a 

greater variation of twist with Mach number than 
did the conventional airfoil. 
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(a) Angle-of -attack (b) Aerodynamic force 
pulse response 

Fig. 3 Transient pulse technique. 


(a) Force coefficients 
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(c) Center of pressure 
Fig. 2 Steady aerodynamic parameters. 



Fig. 4 Transient force for unsteady solution 
with a 0 * 10 held fixed. 64 AMES 
airfoil at M = 0.8. 
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Fig. 12 Flutter results. 


Fig. 9 Unsteady lift Fig. 10 Unsteady 
coefficient moment 

for k = 0.15. coefficient 

for k = 0.15. 



Fig. 11 Flutter root locus for 64 AMES airfoil, 

S = a + 1u), 



Fig. 13 Comparison of flutter boundaries. 
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(a) 64AMES airfoil (b) MBB-A3 airfoil 
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Fig. 14 Effect of elastic axis location on 
static aeroelastlc twist. 
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